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Abstract 

Coupled map lattices of weakly coupled Chebychev maps, so-called chaotic strings, 
have a profound physical meaning in terms of dynamical models of vacuum fluctuations 
in stochastically quantized field theories. Here we present analytic results for the invariant 
density of chaotic strings, as well as for the coupling parameter dependence of given ob- 
servables of the chaotic string such as the vacuum expectation value. A highly nontrivial 
and selfsimilar parameter dependence is found, produced by perturbative and nonpertur- 
bative effects, for which we develop a renormalization group theory. Our analytic results 
are in excellent agreement with numerical simulations of the chaotic dynamics. 
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1 Introduction 



Coupled map lattices (CMLs), as introduced in seminal papers by Kaneko and Kapral some 18 
years ago [TJ [2] , are known to exhibit a rich structure of complex dynamical phenomena [3j HJ [5j 
El [7] . Of particular interest are CMLs that consist of locally chaotic maps. For hyperbolic local 
maps and very small coupling it can be proved [U |9j [TUl [TT] that a smooth invariant density of 
the entire CML exists and that there is ergodic behaviour. However, the case of nonhyperbolic 
local maps, e.g. of local 1-dimensional maps for which the slope is equal to zero at some point, 
is much more complicated from a mathematical point of view, and much less analytic results 
are known [12j [131 CLU EES! EH El EES] • This is the realm of chaotic strings. 

Chaotic strings are 1-dimensional coupled map lattices of diffusively coupled Chebychev 
maps [TH]. They are intrinsically nonhyperbolic, since the local maps have an extremum with 
vanishing slope. They play a very important role in generalized statistical mechanics models of 
vacuum fluctuations [U [TjJ] [20j [2TJ , replacing the noise in stochastically quantized field theories 
by a deterministic chaotic dynamics [221 1231 [21] • These type of chaotic theories are of utmost 
interest in extensions of the standard model where one constructs additional sectors of rapidly 
fluctuating chaotic fields. One possible physical embedding is to associate the vacum energy 
generated by the chaotic strings with dark energy [25] . These novel strongly chaotic models of 
vacuum fluctuations can help to understand reasons why certain mass and coupling parameters 
are realized in nature, others not, and are thus highly important from a fundamental physics 
point of view (see Ref. [4] for more details). 

In this paper, however, we will not deal with the above mentioned physical applications of 
chaotic strings, but merely concentrate onto the mathematics of these nonhyperbolic coupled 
map lattices. There is a highly interesting selfsimilar dependence of the invariant density on 
the coupling constant a of the CML, which can be understood by analytic means and which 
will be the main subject of this paper. 

Let us first recall what has been done so far. In Refs. [IHl [20] weakly coupled iV-th order 
Chebyshev maps [IH], E31 [25], [26], ET] were studied and it was shown that certain scaling properties 
with respect to the coupling a can be calculated perturbatively. However, the fine structure 
of the parameter dependence of important observables of the chaotic string, such as the self 
energy or vacuum expectation value, could not be explained by these perturbative methods. 
In Ref. [28J is was numerically illustrated that the fine structure does occur only for CMLs 
of dimension 1, and not on lattices of higher dimension. In Ref. |29j short chaotic strings, 
consisting only of a small number of lattice points, were analyzed and it was shown that some 
of these short strings (with lengths between 3 and 5 lattice sites) show the same fine structure 
as long strings of the order of 10 000 sites, while for other short strings the fine structure had 
a totally different shape. 

In this paper we report on significant progress to better understand the scaling structure 
of chaotic strings. We will show that the invariant density can not only be understood near 
the edges of the interval on which the local maps are defined, but in the entire interval region. 
We will study the coupling constant dependence of the expectation value of the chaotic field 
variable. This expectation is expressed as a particular integral, which is trivial for the uncoupled 
case (coupling a = 0) but exhibits a highly nontrivial parameter dependence for o ^ 0. There 
are certain invariance properties of observables of the chaotic string under suitable rescaling 
of the coupling that we will describe in terms of a kind of renormalization group theory. For 
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this purpose we introduce two basic scaling functions, which we call 'draft function' and 'blunt 
function', and which describe how parameter changes of the coupling are generating a selfsimilar 
pattern and how certain 'excitations' (travelling maxima) in the invariant density are induced. 
We also further develop the perturbation theory, and illustrate that there are perturbative 
and nonperturbative effects for the chaotic string, with certain analogy to what is known in 
quantum field theories. 

This paper is organized as follows. In Sec. 2 we globalize the perturbative methods developed 
in Refs. [T$| |2"U] and apply them to the vacuum expectation value (0) of the chaotic field 
variable. In Sec. 3 we analyse the integrand leading to (0) and introduce draft and blunt 
functions which can help to analytically describe the self-similar parameter dependence and 
invariance properties of the chatic string. In Sec. 4 we discuss simplified models of the scaling 
structure which are connected to the perturbation theory. In Sec. 5 we look at limiting cases 
for the partition integrals that occur when the vacuum expectation is calculated. In Sec. 6 we 
provide analytical expressions for the (exponential-type) draft function. Finally, in Sec. 7 we 
draw our conclusions. More details on explicit calculations can be found in the Appendices. 



2 The perturbative method 

As shown in Ref. [20] in detail, the scaling behaviour of the self energy (and interaction energy) 
of the chaotic string based on Chebyshev maps T N ((p) of order N is directly related to the 
scaling behaviour of the distribution function (the 1-point invariant density) p(0) of the string. 
For the chaotic string of type 2B (cf. Ref. [I]) with iterative prescription 

C +1 = (i - a)T 2 (C) + a -{^ +1 + <j>U) =■■ T 2a (0r; 0? +1 + <t>U) (i) 

the distribution function close to the upper boundary 0=1 turned out to be a sum p a {l — ax) w 
E P P =T Pi p) (l - an) of functions p& (in the following sometimes called p-iterates) with0 

- ax) = , 1 f tott+W+Pott-W- (2) 



7iy2a(l - a) J Ap^Jx/Ap - rf(0 + ) - r 
while close to the lower boundary = —1 we obtain p a (ay — 1) ~ p£'(ay — 1) with 

p(P)( y j) - 1 f PO^+^+PO^-)^- (3) 



■Kpa{l-a)J v /y-(0 + + 0_)/2 

Po(0) — l/7T\/i — (j)' 2 is the invariant density for the uncoupled string (a = 0). We use the 
notation + := 0_ := 4> '■= 0™, and 

W) := o2. — ^ — • ( 4 ) 

* q=0 JV 



1 Compared to Ref. [20] we use a different definition (j4]) for the function r p N (<j)) in order to keep it positive. 
Note that by convention the integration range is given by the positivity condition for the radicand 20 . 
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Both approximations can be formally combined to the formula 

Pmax 

pM= Ep1 p) (0) (5) 

p=0 

This approximation is valid only close to the boundaries = ±1. This is shown in Fig. [T] (close 
to the lower boundary) and Fig. [5] (close to the upper boundary). For the main part of the 
interval G [—1,1], however, there are large deviations between numerical results obtained by 
iteration of the CML and the above formula, regardless of the order p max of the approximation. 
Nevertheless, a main result of this paper is an extension to a formula that is valid for general 
values of 0. 



2.1 Extension to the whole range 

The key observation for this extension is given by the insight that the integrand of 

p(o)(0) = 1 f p ((j) + )d(f) + po((j)_)d(j)_ ^ 



nVl-aJ ^/2(0 + l)-a(0 + + </»_ + 2) 

is given by the derivative of the inverse map of the map (JTJ) 

1 d 



2^/T^a~J2{<j) + 1) - a(0+ + 0_ + 2) 



7^(0; + + 0_) =: T^Ofc + + 0_). (7) 



Of course, the inverse T 2 ~ x (0; 0+ + 0_) of the function T 2a (0; 0+ + 0_) given by the quadratic 
equation (JTJ with T 2 (0) = 20 2 — 1 is not unique. We have used the positive square root 
T 2 +(0; + + 0_) according to the notation 



T 2 t(0;0 + + 0-) := ±« 
in order to write 



pfW = 1 I Po(0+)^ + Po(*-)#-^7r+(^; </>+ + </>-). (9) 



7T 

Proceeding to the first iterate, we extend to 

2 r 

Pa } (0) = / Po(0+)d0+Po(0-)d0- X 

7T J 

x (T 2 -(0; T 2 (0 + ) + T 2 (0_)); + + 0_) . (10) 

For details of this and the following see Appendix A. For an arbitrary order p we finally arrive 
at 

Pi P) (0) = / Po(0+)d0 + Po(0-M0- X (11) 

7T J 

x -4t 2 + {T 2 + a (• • • (0; T 2P (0 + ) + T 2P (0_)) . . . ; T 2 (0 + ) + T 2 (0_)) ; + + 0_) . 
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Figure 1: Numerical result obtained by (a) long-term iteration of the CML and (b) the analytic 
approximation (jSJ) at maximal order p max = 7. The figure shows the rescaled distribution 
function (p a ((f>) — po (</>)) I Pq{4>) a ^ coupling a = 0.00007625 close to the lower boundary = — 1 
in dependence on log 4 (l + </>). 
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-log 4 (1 

Figure 2: Same as Fig. 1 but close to the upper boundary = +1 in dependence on — log 4 (l- 
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2.2 The "path of roots" 

The path "+ ■ • • H — " for selecting the nested square roots turns out to lead to an excellent 
approximation close to the upper limit <fi = +1. Of course we have to ask why this path needs 
to be taken and not a different one. The answer to this question is related to the trajectory 
that a starting value close to zero takes under nearly unperturbed iterations. Starting from 
= we obtain T 2 (0) = -1, T 2 (T 2 (0)) = +1, T 2 (T 2 (T 2 (0))) = +1 and so forth. This is just the 
time-reverse of the relevant "path of roots" . However, if we want to extend the approximation 
of the invariant density in the way mentioned above, we need to take into account all paths. It 
turns out that a great amount of other paths returns contributions of lower order. Moreover, 
these contributions add up in a way that a general factor 2 /it is factored out. The final result 
(understood as sum over all paths of roots) 



p O (0+)d0+Po(0-)#- X 



(12) 



d 



x — T," 1 T,- 1 fc 1 (. . . ; T 2 (T 2 (0+)) + T 2 (T 2 (0_))) ; T 2 (0 + ) + T 2 



+ 



for a single p-iterate (without adding lower iterates) approximates the distribution function to 
an increasing degree with increasing order p. 
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Figure 3: Comparison of the rescaled p^((fi) (p = 5,6,7) with the numerically obtained his- 
togram result. In this plot the histogram is shifted by the amount +0.2 for better visibility. 



In Fig. [3] the coincidence of higher p-iterates with the distribution function is shown. For 
the coupling we have chosen a rather high value a = 0.00122 as we are no longer interested in 
the perturbative region close to the boundary but in nonperturbative effects in the middle of 
the interval. As compared to Fig. 1 and El the position and magnitude of the local maxima 
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(often called 'excitations' in the following) are mirrored in a precise way. Even the excitation 
close to = —0.15, i.e in a highly nonperturbative region, is exactly modeled. 

The excitation at ~ —0.15 is special. Indeed, if we look at a sequence of diagrams similar 
to Fig. 3 where the coupling a gradually decreases, then this excitation moves to the left while 
all others move to the right. These and other counter movements account for the so-called 'fine 
structure' of observables of the chaotic string. By this we mean a nontrivial local a-dependence 
of observables as illustrated e.g. in Fig. 4 and 5. The reason for the different behaviour of the 
above excitation is that this excitation belongs to a different "path of roots" than the ordinary 
one ("+■•• H — "). It can be associated with a process that is mirrored at the lower boundary. 
This becomes apparent when we look at a sequence of diagrams as in Fig. 3 again. At the 
moment when the extraordinary excitation disappears at = —1, a new ordinary excitation 
appears at this point and moves to the right. 

2.3 Expectations of observables 

The fine structure we are trying to understand in this paper is the a-dependence of observables 
(functions of 0) such as e.g. the self energy V(</>) = — 20 3 /2 calculated with respect to the 
distribution function p a {4>) (cf. Fig.H]). These types of expectations of observables as a function 
of a turn out to be of utmost interest for the physical applications described in Ref. [I] . It turns 
out that the vacuum expectation value (0) of the string leads to the same (but less pronounced) 
fine structure, as it is shown in Fig. [5j Therefore, in the following we will concentrate onto the 
calculation of (0) = J 0p a (0). Approximating p a {4>) by the p- iterate p^(0), one can perform a 
{p + l)-fold integration by parts to obtain 

(0) a = J 0pi p) (0)# = / Po(0+)rf0+Po(0-)rf0-T 2a (T 2a (T 2a ( . . . T 2a (0^; + + 0_) . . . ; 
T 2P - 2 (0 + ) +T 2P - 2 (0_));T 2P -i(0 + ) +T 2P -i(0_));T 2P (0 + ) + T 2P (0_))d0^. (13) 

This equation reminds us of results obtained in Ref. [29] for an open string of length 3@ In 
order to better understand how the excitation can switch to a different "path of roots" if a is 
changed we take a detailed look at the integrand in Eq. fflBl as a function of ( - p - ) (renamed as 
0) for different values of (0 + ,0_). For (0 + ,0_) = (+1, +1) the integrand at a = 0.00122 is 
shown in Fig. [HI We see that while the minima of this strongly oscillating function take the 
value —1, some of the maxima take values less than 1. The most extreme situation is found 
at = where the maximum is on the way to decay completely if a changed in a significant 
way. When this happens, the two neighbouring minima unite into one new minimum. From the 
calculations in Refs. [T9"l [SO] it can be seen that the height of maxima is related to excitations 
of the distribution function. The decay of a maximum at = for increasing values of a 
thus means that an excitation ceases to exist at the lower boundary. The vanishing of the 
extraordinary excitation mentioned before is related to such a formation of a united minimum. 

2 Note that due to a numerical error, in Ref. [29] is was erroneously stated that the open string of length 3 
does not reproduce the fine structure. This will be corrected in an erratum to Ref. [29] . 
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1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 

a [10" 4 ] 

Figure 4: Expectation value of the self energy V(<f>) in the range a G [0.00015, 0.00060] for a 
2B-string of 10 000 lattice points and 10 000 iterations 
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Figure 5: Expectation value of <p in the range a e [0.00015,0.00060] for a 2B-string of 10 000 
lattice points and 10 000 iterations 
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Figure 6: Integrand of Eq. ({TBI) for p = 5 at a = 0.00122 and 0_ = <p + = 1 

3 Structural changes of the generating partition 

The change described before is a change of the topological structure because one of the local 
oscillations which partitions the interval [—1,1] into TP intervals vanishes. Indeed, the different 
partitions stand in direct correspondence to the different roots of the integrand. For vanishing 
coupling a = the maxima and minima are located exactly at sin(7rt/2) for the values t = 
±2n/2 p and t = ±(2n± 1)/2 P , resp., where n = 0, 1, . . . , 2 P_1 counts the extrema. The minima 
divide the interval [—1,1] into TP intervals which we call the elements of the generating partition 
(see Fig. 6). If one of these elements vanishes, the topological structure of the partition changes. 

One can use the generating partition to split up the integral in Eq. ffl3|) into a sum of TP 
partition integrals. From this point of view, the decay of a maximum means the vanishing of 
the corresponding partition integral. In the following we take a closer look on how an increasing 
coupling a influences the contribution of the integrand through the different elements of the 
generating partition. In order to distribute the elements of the generating partition uniformly 
onto the interval [—1,1], we use a new variable t, with = sin(7rt/2). 

3.1 Decay of the maxima 

A general feature of the integrand is that maxima which for a = are located at equally spaced 
positions ("starting positions") to decay in a similar way when a is increased. This property 
forms the basis for our renormalization group approach in the following. Let us assign a class 
k to certain subsets of maxima as follows: The maxima with starting positions t = and 
±1 are of class k = 0, the maxima with starting positions to = ±1/2 are of class k = 1, the 
maxima with starting positions to = ±1/4, ±3/4 are of class k = 2, and so on. In general, the 
starting positions of class k are given by to = ±(2n + l)/2 fc , n = 0, 1, . . . , 2 fc_1 — 1. We note 
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the similarity with the definition of /c-cylinders in the thermodynamic formalism of dynamical 
systems [30J. For (0 + , </>_) = (+1, +1) the height of a maximum of class k is given by 

b = cos( v / ^2 p - fe+1 ). (14) 

From Eq. (fl4|) we read off that for a specified scaling interval where a increases by a factor 
4 (e.g. a e [0.00015,0.00060]), precisely one single specified class of maxima (and, therefore, 
partition integrals) disappears, and this disappearence takes place synchronously. Apparently 
the successive decay of maxima produces global scaling features, such as the invariance under 
the transformation a — > 4a, but it is not responsible for the fine structure within a given scaling 
region, which is more related to infinitesimal changes of a. To better understand the fine 
structure, in the following we will introduce two important auxiliary functions: Draft functions 
and blunt functions. 

3.2 The draft function 

If elements of the generating partition vanish, the neighbouring partition elements fill the 
gap, being drawn close to the position where the topological change happens. The difference 
between the starting position of a given maximum and the position for a ^ turns out to be 
a linear function in a. The slope, therefore, is a well-defined function called draft function, 
d(4>o) = d(p a (a; (p )/da where </> a (a; </> ) is the position of the maximum with a (O; 4> ) = (f> . For 
large p and k this is a fractal function, and the basic features visible in a plot are independent 
of p and k. The function is shown in Figs. U\ E]and[9]for (0 + , <p_) = (+1, +1), (+1, —1) (which 
is the same as for (—1,-1-1)) and (—1,-1). We found an analytical expression for the draft 
function. With O — sin(7rt /2) the analytical formula for the draft function of class k reads 

dWfti; 0-) = \ {d (k) (t ; <j> + ) + J^(t ; <P-)) (15) 

(the tilde indicates the linearization of the first argument by replacing <fr = sin(7rt /2) by t ) 
where 

fc— l 

d {k) (t ;<p±) = o^(to) + Y,TA-<P±)d < f\t ) (16) 

1=0 

and 

S k \t Q ) = 2-*d?}(to) - E (^- k+l + cos(2'7rt Q )) d) k \t ), d^ito) = . 1 (17) 

The draft function d {k \t ; <f>+, 0_) is smooth for t = ±(2n + l)/2 fc , n = 0, 1, . . . , 2 k ' 1 - 1 but 
singular for values to corresponding to a class smaller than k. 

3.3 The blunt function 

Eq. (|14|) is valid only for the specific value (<f) + , </>_) = (+1, +1) of the neighbours. For general 
neighbouring values <p + , 0_ the height of maxima of class k is given by 

b p , k (<P + , <P-) = cos (V- fe+ Varf (-T 2 *(0+), -T 2 *(0_))) (18) 
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Figure 7: Draft function for 0+ = 0_ = +1 
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Figure 8: Draft function for <p + = +1 and 0_ = — 1 
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where r p N is defined in Eq. PJ. We call b Pt k the 'blunt function'. Both draft and blunt functions 
are important for our renormalization group description in the following. Basically the draft 
function describes how positions of maxima change under parameter changes, and the blunt 
function describes how the values of the maxima themselves change. 

If the integrand for different values of (</>+,</>_) is affected by a relative draft, the blunting 
is accelerated because after integration over p o (0+)po(0-) the integrand (in <f>) is smeared out, 
giving rise to a premature breakdown of the integral. 

4 Renormalisation group approach 

The simultaneous parameter-dependent deformations described by draft and blunt functions 
are the basis for our treatment in the following. To understand the complex a-dependence of 
observables such as the vacuum expectation value, we need to evaluate the integral ( [TBI) . This 
integral, as mentioned before, is split into many partition integrals. These partition integrals are 
defined as the part of the integral ({TBI) where the integration range is restricted to the interval 
given by the positions t a and tj, of two adjacent minima surrounding a decaying maximum with 
starting position t , 

(4>)M = ~£f(t)<™$)dt (19) 

where the integrand is given by 

f(t) := J p o (0+)#+Po(0-)#-T 2a (. . . T 2a (sin(7rf/2); 0+ + 0.) . . . ; T 2P (0 + ) + T 2P (0_)). (20) 

In Figs. [10] and [TTJ we show a sequence of partition integrals and their sum for values of p and 
k with constant difference p — k = 6. At first sight one might conjecture that all these partition 
integrals vanish at the same maximal value of a, called the maximal reach for a given (p,k). 
However, our careful analysis shows that due to different values of the draft function, many 
partition integrals vanish at values below the maximal reach for a given (p, k). 

Even though the draft function does not depend on p and k, the effect of partition integrals 
to vanish below the maximal reach increases with the order because the maxima are closer 
together and are smeared out faster. For this reason, also the maximal reach itself is decreasing 
in the sequence of the two diagrams shown in Figs. [10] and [TTJ 

In the spirit of a suitable renormalization group theory, we are now looking for a transition 
p — > oo that preserves the shape and sum of the partition integrals. Comparing the case p = 14, 
k = 6 as displayed in Fig. [T2] with Fig. [TU] (p = 10, k = 4), and the case p = 13, k = 6 in 
Fig. EH with Fig. [TTJ (p = 11, k = 5), we find that the quantity to be kept constant is given 
by the reduced order p = p — 2k which in the two examples is given by p = 2 and p = 1, 
respectively. However, in Figs. [12] and [13] we have stretched the axis for a by a factor of 4 
and 16, respectively. Because we are looking at the behaviour in the scaling region only, this 
means that instead of the coupling a one should use the reduced coupling a = 4 p ~ k a which for 
the intervals shown in Figs. [TUHT21 is found in the interval a e [0,4]. With the above rescaling 
prescriptions, the overall problem is formulated in a scale invariant way in the limit p — > oo. 
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Figure 10: Contributions and sum of all the partition integrals for p — 10 and k = 4 that vanish 
in the displayed interval for a. For better visibility, the sum is multiplied by — 2 1 ~ fe . 
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Figure 12: Contributions and sum of all the partition integrals for p — 14 and k = 6 that vanish 
in the displayed interval. For better visibility, the sum is multiplied by — 2 1 ~ k . 
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Figure 13: Same as Fig. 12 but for p — 13 and k = Q. 
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4.1 The model function 



The Feigenbaum fixed point function in Feigenbaum's renormalization group approach for 1- 
dimensional critical systems can be well approximated close to the vicinity of the maximum of 
the map [301 Ell E2]. Similarly, we can also find a good approximation of the integrand in our 
problem if we are close to a decaying maximum. We call this approximating function of the 
integrand the 'model function'. Using the analytical expressions for draft and blunt functions, 
the integrand in Eq. (|2"U|) in the region close to a decaying maximum is approximated by the 
model function 

/(At; to) = y p o (0+)rf0+Po(0-)#- x 

x cos (^vr 2 (At - arf»(t o; cf> + , 0_)) 2 + 4 1 ~ fc ar 2 (—T 2 * ; (</>+), -T 2fc (0_))^ (21) 

which is derived from Eq. (1TB]) . This turns out to be a quite precise approximation for f(t + At) 
in Eq. (|2"0"]) where At is the deviation from the starting position t of the specified maximum. 
If in addition to the reduced order and coupling one uses the reduced deviation At = 2 p At e 
[—1,1], the limit p — > oo can be performed for the model function ( I2ip . rewritten as 

f^(Ai;t ):=f(2- p Ai;t ) = [ p (0+) dp +Po (0_) X 



x cos ^tt 2 (At - 2Parf»(t ; <P+, <j>~)) 2 + Tarf (-T 2 ,(0 + ), -T 2 ,(0_))j . (22) 

The function r 2 5O (0 + , 0_) uses the fc-fold Chebyshev map for the arguments <p± while the argu- 
ments of the draft function are just (f>±. Keeping the reduced order p = p — 2k constant, the 
class k increases together with the order p. Because of this, the high-frequency blunt can be 
separated from the low-frequency draft, leading to an effective blunt of 8/3 (see Appendix B). 
The final result for the model function in the limit p — » oo reads 



/<«(At; t ) = / M0+)#+Po(0-)#- cos ^vr 2 (At - 2^arf»(t ; + , 0.)) 2 + . (23) 
4.2 Fourier transformation 

We note that the integrand in the above equation ( [23]) is essentially of the form g(a,t) : = 
cos(\/t 2 + a). Inspired by quantum field theory, we proceed to a Fourier transform denoted as 

/oo 
cos(vt 2 + a) cos(ut)dt. (24) 
-oo 

It can be shown that g(a,u) = for u 2 > 1, but otherwise 



g(a, u) = 2nS(l - cu 2 ) - ™ J x Ua{\~uj 2 )] , 

^0(1 -u 2 ) ^ / 



(25) 
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where J\(z) is the Bessel function of the first kind. One obtains 

a du 



g(a,t) = cost — / — = =Ji ( \/a(l — to 2 )] cos(ut). (26) 
J - 1 2^a(l-u 2 ) V v J 

If this is inserted into Eq. f|23|) . the model function can be expressed as 



f®{At;to) = J p o (0 + )#+Po(0-)#- 
4a du 



cos 



(i-K (At- 2^ad (fc) (t o ;0+,0_ 



Ji (^8a(l -u; 2 )/3J cos (mw (At - 2^aa T(fc) (t ; 0+, <£_)) 



] I l/OIII I. I I'llh I 2l V — _ ' llll I Ifl- <.'-'-■ <■>- IJ I . (27) 

3^/8a(l - w 2 )/3 ' 
4.3 Bessel integrals and Feynman webs 

The occurence of the Bessel function of degree one discloses a surprising connection to the Feyn- 
man diagrams of elementary particle physics. In calculating Feynman diagrams in configuration 
space and afterwards returning to momentum space with Euclidean metrics in D dimensions, 
the exponential function can be expanded into Gegenbauer polynomials Cf(w) [33J, 

ip-x ~nf \ \ I r I ;lf\ i /\ j ( \/^yA / P x 



e- = r(A)(^) E^A + O^WCM^ , (28) 

where D = 2 A + 2. Note that p ■ x = p^x^ is the scalar product of two D-dimensional vectors, 
while px is the product of the two absolute values p = y/p^ and x = ^/x^x^. If the integrand 
does not depend on the D-dimensional vectors but only on their absolute values, then the 
D- dimensional integration measure of the Fourier transformation simplifies to 

J e ip x d D x = 2tt a+1 X J x (px)x 2X+1 dx. (29) 

An interesting observation is that the model function contains the Bessel function (and a similar 
factor in front of it) of degree A = 1, which corresponds to-D = 2 + 2 = 4 space-time dimensions. 
If one aims for a physical embedding of the chaotic string dynamics into generalized versions 
of quantum field theory [I], then this means that A = 1 cites four- dimensional space-time 
without the explicit emergence of this space-time itself . If we interpret the integrand of the 
Bessel function as a correlator function in (obviously, four-dimensional) configuration space, 
the integration over <p + and <ft- basically means the insertion of two vertices into the diagram 
which are connected by a propagator, i.e. a virtual particle transition. The fish-type diagram 
in Fig. [HI however, is a compactified version of the Feynman web that was used in Ref. jl] to 
give a physical interpretation to the chaotic string dynamics. 



5 Limiting cases 

Clearly the exact renormalization group treatment requires the limit p — > oo. On the other 
hand, of relevance for a scale invariant formulation is the reduced order p = p — 2k. In the 
following we look at different limiting behavior depending on the relative size of p and k. It 
turns out that one can distinguish two different limit ranges which we call the perturbative and 
the nonperturbative range. 
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5.1 The perturbative range 

Increasing the reduced order p, the maximal reach does not increase in the same way. Instead, 
more and more partition integrals vanish close to the maximal reach (see Fig [15]) . This means 
that compared to the blunt effect, the draft effect can be neglected. Note the similarity of 
the three displayed pictures in Fig. 15, which is due to the fact that the reduced coupling 
a = 4 p ~ k a is used. This range of parameters p and k does not contribute to the fine structure 
but is responsible for the global selfsimilar behaviour of expectations of observables. 

5.2 The nonperturbative range 

More interesting for the fine structure is the case where the reduced order p becomes negative. 
Fig. [T6l unveils a rich structure. While Fig. [16] is calculated by using the model function ( 1271) . 
the calculation for the open string of length 3 leads to the same result. Finally, the result for 
class k = 9 (with the same reduced order p = —3, or p = 15) is again very similar. This means 
that the model function correctly describes the effects seen in this nonperturbative range. Very 
characteristic is the kink close to a = 0.5 which is relevant for the fine structure. This kink has 
its origin in a majority of partition integrals swinging twice before they vanish. 

5.3 Development of bounded states 

Characteristic for the nonperturbative range is that for a given decaying maximum, the neigh- 
bouring maxima are influenced strongly by the decay of the central maximum. They decay 
nearly as fast as the central maximum. If we use as an analogy the language of particle physics, 
then maxima seem to interact and to formate something like a bounded state. An interesting 
open question is whether these "bounded states" can be described properly within our current 
formalism based on partitions, or whether a different approach is necessary for this. 
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Figure 15: Contributions and sum of partition integrals of class k = 7 for increasing order 
p = 17, 18, 19 (from top to bottom, corresponding to the reduced order p = 3,4, 5) as functions 
of a. For better visibility, the sum is multiplied by — 2 1 ~ fe . 
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Figure 16: Contributions and sum of partition integrals for p = 13 and k = 8 as a function of 
a G [0,4]. For better visibility, the sum is multiplied by — 2 1_fc . 



6 The exponential draft function 

In the following we will perform further calculations on the draft function, which will ultimately 
lead to the implementation of graph theoretical methods previously introduced in Refs. [341135] . 
Let us define an exponential draft function by 

h{u- 1 ) := J* exp (-icuS k \t ; 0+)) p (<P+)d<P + . (30) 

With the help of this the model function can be written as 
/®(At;t ) = Re{e i7rAi ~h 2 {2- i '- 1 7Ta-t 



o) 

[ +1 ; AadU Ji (jSa(l-uj 2 )/3) Re (e i ^ Ai h 2 {2^- 1 Tiauj ] t j) . (31) 
•/-i 3^(1 - w 2 )/3 V V 7 



The relative minima (and maxima) of the model function which define the boundaries of the 
partition integrals are found by calculating the first derivative with respect to At and deter- 
mining the zeros of this derivative. One obtains 



Im (e^ A *l 2 (2^- 1 7ra;t ) 

= t ; 8hudU — jJjSajl-u^^lmU^ Ai h 2 (2~ ii - 1 7rdu J ;to)) , (32) 
3^(1 -u 2 ) /3 V v J v ' 

using the fact that the integrand is an even function of uj. 

The exponential draft function can be visualised by a parametric curve in the complex 
plane, as it is shown for k = 8, p = 13 and a particular partition integral specified by n = 42 
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Figure 17: Parametric curve of the exponential draft function h(2~P~ 1 nauj; t ) for k = 8, p = 13 
and n = 42 with reduced coupling a = 2.4 and u> G [0, 1]. The twelve points represent 12 
equidistant values of ui for a fixed value of a — 2.4. 
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in Fig. [T71 The long range of the curve is due to the fact that 2~^~ 1 7ra = And is not a small 
number. Note that the curve for u> G [—1,0] is obtained by mirroring at the real axis. The 
parametric curve shows quite an irregular behaviour which is hard to predict. 



6.1 Analytical expression for the exponential draft function 

One can perform a concrete calculation to obtain the exponential draft function using graph 
theory. Inserting Eq. ( Tl6|) into Eq. ( )30|) . one obtains^ 



h(u);t ) 



k-1 



exp ( -iuJ2 T 2i(<P+)d ( l k \t ) ) p(<p+)d<p 



1=0 



00 I A, ,\n d k— 1 



E 

n=0 



111 



^r 2i (0 + )rf? ) (to))'po(0 + )rf0 + . 



(33) 



1=0 



In the following we use the short hand notations d := d^ k \to) and di := d\ (to) and the 
convention that sums run over values where all indices under the sum take nonnegative values 
less than k. By performing the integration over + one obtains 



h(u;t ) 



+ 
+ 



+ 



+ 



-iu) A ^ 2 Mo;) 4 

2 , «J «2+l«Z+2 + 



16 



32 



(~^) z 
64 



£4 



32 

— io;1 



E dqdi + id[ +2 di +3 + 



-IUJ 



E <tfdf +1 



(-iuj) 



32 £#fc + i£4 + ^£#*. 



Z+2 



— zo; 



48 ^ 1 + 64 



64 

-iw) 6 ^ 2 3 (-ZU;) 6 ^ 4 (-ZW) 6 ^ 3 

2^ «i+i««+3 H 2^ "i ««+2««+3 2^ "i d l+l d l+2 



384 



96 
-iuj 



dfdi+id. 



1+2 



768 

Ho;) 6 

128 



^c/^ +1 rff +2 -ll- 



-20; 



1536 



E d i d i+i 



k-l 



IUJ 



64 128 



£W) +^E^ 



256 



384 



E^ 2 + ■■■ 



(34) 



Expressing the higher-order correlation of (uncoupled) Chebyshev maps in terms of a graph 
theoretical approach [3U [35] , the various contributions above represent possible integer parti- 
tions of powers of 2 and can be represented as binary forests. It can be shown that the series 
can be resummed as an exponential function h(cu;t ) = exp(ih z (cu;t )) where the exponent 
h z (uj;to) contains only trees, 

r / x • , ( - i^) 2 ^ , 9 (— iu) 3 ^ , 9 , 



3 Without loss of generality, T 2 i(—(f>±) can be replaced by T 2 i(4>±) under the integration. 
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Ho/) 4 v 2 (-io;) 4 v 4 (-w) 5 v 2 

Hwf v 4 (-zu;) 5 v 2 3 (-w) 5 v 4 

(-iw) 6 ^ 2 (-iw) 6 ^ 2 3 ( - ^) 6 V^j4j 
H g^ d i «/+l«/+2««+3««+4 H 2^ d / «Z+l«Z+3 H ^gg 2^ d « «Z+2«Z+3 

(-^) 6 V , 2 ,3 ^ (-io;) 6 v 4 Hfv^ r/ 3 

2^ «Z a i+i a l+2 777 2^ a i «Z+l«Z+2 77777 2^ a i a l+l a i+2 



4 g ^ i i+i ^ gg ^ i *tl 128 

1536 ^ ' m 576 ^ ' V ; 

6.2 Single trees and multiple trees 

The idea of using binary forest graphs for our calculation can be summarized as follows: Let di 
be the coefficient of the Chebyshev polynomial T 2 ;(0 + ). Under the replacement <p + = sin(7rt + /2) 
one has T 2 o(0 + ) = sin(7rt + /2), T 2 i(0 + ) = — cos(7rt + ), and T 2 ;(0 + ) = cos(2' _1 7r£ + ) for I > 1. 
Therefore, under this replacement, the Chebyshev poynomials can be represented using complex 
exponential functions exp(±2'~ 1 i7rt + ). In calculating the integral in Eq. ( 13"3"j) . a product of 
exponential functions corresponding to the coefficients di +i in the sums in Eq. is integrated 
over cf) + . However, the integration of this product gives a non- vanishing contribution only if 
the exponent of the complex exponential vanishes (for details, cf. Refs. [33] 135]). 

To give an example, the /-th contribution to dfdi+i in Eq. is given by integrals of 
the exponential function 

exp ((±2' ±2 l ± 2 l+l )mt + ) = exp (2'(±2° ± 2° ± 2 l )mt + ) . (36) 



This integral gives a non-vanishing contribution only in the case where ±2° ± 2° ± 2 1 = (all 
the signs being independent). The only possibilities in this case are 2° + 2° — 2 1 = 1 + 1— 2 = 
and —2° — 2° + 2 1 = — 1 — 1 + 2 = which are different only by a global sign. One can write 
2° + 2° = 1 + 1 = 2 = 2 1 which is a unique partition of the highest power. In general, the 
binary sum condition is given by 

i n m 

E E ^ m = (37) 

m=0 n=l 

where s n = ±1 are the possible signs in the binary sum condition. The same powers of 2 in this 
condition are represented by points (leaves) in the same layer of a binary forest of graphs. In 
order to understand these graphs and their construction methods, one starts with the highest 
power. The highest power 2 l can be decomposed into lower powers by using 2* = 2 1 " 1 + 2* _1 . 
This is represented by two lines going down from this (highest) point of the graph. If both lower 
powers are found in the binary sum condition (like in 2 1 = 2° + 2° in the example above), the 
lines are closed by points (leaves are drawn) and a tree is created. If one (or both) of the lower 
powers are not found in the vanishing sum, this lower power again is decomposed in smaller 
powers, until the powers are met by contributions in the binary sum condition. 

The main part of the exponent h z (u]t ) contains single trees, corresponding to unique 
decompositions of the highest power of 2. They are given by 

n = 2 : V] df means 1 — 1 = • 
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n 



n 



3 : ^2 dfdi + i means 1 + 1 — 2 = « // *\ 

4 : dfdi +1 di +2 means 1 + 1 + 2- 4 = 



n = 6 : 



n = 5 : dfdi +1 d t+2 di +3 means 1 + 1 + 2 + 4 — 8 = 



^2 dfdi + 2 means 1 + 1 + 1 + 1- 4 = 





E j2 



dfdi + idi +2 di +3 di + 4 means 1 + 1 + 2 + 4 + 8 — 16 = 




^ dfdi +2 d i+3 means 1 + 1 + 1 + 1+ 4- 8 = 



^ dfdf +1 d l+3 means 1 + 1 + 2 + 2 + 2- 8 = 

The coefficients for these contributions are given by 

1 





)n— 1 



n !ni! • • -n n _ 2 ! 



(38) 



(39) 



where nk counts the occurence of 2 k in the decomposition. However, there are also multiple 
trees, like for instance J2dfdf +1 , corresponding to +1 + 1 — 2| + 2 — 2 = where the vertical 
bar separates the two different independent relations +1 + 1 — 2 = and +2 — 2 = 0. In order 
to deal with these contributions and their coefficients, it is appropriate to work with so-called 
summation ordering (for more details see Appendix C). 



6.3 Summation ordering 

If for multiple sums one excludes repeating indices, one ends up with a summation ordered 
expression, similar to time ordered expressions in quantum field theory. One obtains 



16 



+ ^ £ d] ± 4 + J24 + 4^ E dld l+ldl+2 d l+3 



+ 
+ 



32 



-ltd) 



32 



32 
96 



E tfdi+i ^2 $ 



-10J) 



384 



E rf M+2 + 



-IU! 



64 



E dfdf +1 



4 (—iio) 6 ^ 2 (—iio) 6 ^ 2 3 

2^ d t d i+ i H — — 2^ d l di +1 d i+2 di +3 di +4 H 2^ d l+1 d [+3 



384 
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-iu) G ^ 4 (-zo;) 6 ^ 2 3 (-iw) 6 ^ 4 

2^ «/ ai +2 cti +3 H — — 2^ "/ d i+i a i+2 H 77^— 2^ "z «i+ia/+2 



' 768 96 ^-i-w'-r- 192 

H jt^t- 2^ "z «J+i«z+2 + 1536 "i "z+i H — Z^ "J «z+i«z+2 2^ % 

+ ^ E dld l+1 ± 4d P+1 + { -^f £ d? 

+ ^E^E4 + ^E^E4E4 + ---) = (40) 

= 5 exp ^ - + £ d] + £ d 2 d m + £ d?d {+1 d, +2 

(-iw) 5 ^ 2 3 (-iw) 5 ^ 4 (-iw) 6 ^ 2 7,, 

H "Z "Z+l H g~g~ Z^ "Z "'+1 H Z^ "Z «Z+l«Z+2«Z+3«Z+4 

1 H^V ^b 1 (~^) 6 V 747 , 1 (-^) 6 y^2, 3 
H Tj^r— Z^ "z "z+i«z+3 H Z^ "z «i+2««+3 H ^ — 2^ "z «;+i"i+2 

1 H")V ^ H 1 H^) 6 V . 2 7 ,3 , H^l!v^ 2 

H Z^ "z «z+i«z+2 H 2^ "z «z+i«z+2 + 1536 2^ "/ "2+1 

2304 ^ ' y v ; 

Here a prime at the summation symbol means disjunct indices, the symbol S indicates the 
summation ordering. Summation ordering defined here, therefore, means that in a product of 
sums the indices are disjunct. As the time ordering cannot be used directly for the exponent 
of a time ordered product in quantum field theory, the same holds for the summation ordering: 
it is defined only if the exponential series in Eq. f j4TT) is expanded as in Eq. ( l4*0|) . However, 
the coefficients are all positive, and by combinatorical means (worked out in Appendix C), the 
coefficients of the non-ordered series can be retained. 



7 Conclusion 

In this paper we have presented novel analytic results for an important class of coupled map 
lattices (CMLs), so-called chaotic strings, which play an important role in extensions of stochas- 
tically quantized field theories. We were able to derive an analytic expression for the invariant 
1-point density of the coupled map lattice as a function of the coupling, extending previous 
results obtained in Refs. [TP] [2D] . The invariant density was approximated to p-th order in an 
iterative recurrence scheme, giving exact result for p — > 00. 

This scheme allowed us to better understand the nontrivial dependence of expectations 
of observables (such as the vacuum expectation value) on the coupling constant. A complex 
selfsimilar structure was found, for which we developed a renormalization group theory. Two 
types of scaling functions were introduced (draft functions and blunt functions) which allowed 
us to comprehensively describe topological changes, as well as gradual changes, of the generating 
partition as a function of the coupling a, thus understanding the corresponding changes of the 
integrals that occur when expectation values are calculated. 
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Analytical expressions were found for both draft and the blunt functions. Based on these, in 
the vicinity of local maxima the scale and parameter invariance could be described in an analytic 
way by introducing a renormalization group 'model' function (the analogue of Feigenbaum's 
fixed point function for 1-dimensional critical maps). By implementing a graph theoretical 
method, explicit calculations could be performed for this model function, with surprising con- 
nections to Feynman graphs in quantum field theory and other graphs (binary forests) developed 
to understand the higher-order correlations of uncoupled Chebychev maps. The exponential 
draft function as the main element of the model function was expressed as a (summation or- 
dered) exponential series. 

Our results significantly advance the analytic understanding of chaotic nonhyperbolic CMLs, 
of relevance in many different areas of physics, and are in excellent agreement with direct 
numerical simulations of the CML. 
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A Explicit calculation for the first and second iterate 

In order to explain our method to obtain the first iterate in Eq. (jTUl) . we calculate 
d 



Tt (T+ (0; T 2 (0+) + T 2 (0_)) ; + + 

= T£ (0; T 2 (0 + ) + T 2 (0_)) T£ (T+ (0; T 2 (0 + ) + T 2 (0_)) ; + + 0_) . (Al) 
Close to = 1 we obtain 

T+'(l-ax;T 2 (0 + )+T 2 (0_)) = ~ + 0(a), 

T+(l-ax;T 2 (0 + )+T 2 (0_)) = 1 - ax' + 0(a 2 ) (A2) 

where 

*' = ^ + i(T 2 (0 + )+T 2 (0_)-2). (A3) 

Inserting the positive square root T 2 +(0'; + + 0_) « 1 — ax' for 0' into T 2 ~ 1 '(0'; + + 0_) does not 
lead to the approximation (J2J). However, if we insert the negative square root T 2 ~(0'; + + 0_) ~ 
— 1 + ax', we obtain 

r+'(-l + ax'; + + 0_) » 1 (A4) 

2y/2a(l-a)Jx/4 - ^(0+) - ^(0-) 

and together with T 2 ~'(l — ax; T 2 (0 + ) +T 2 (0_)) f» —1/4 the integrand of the first iterate, which 
is now given by 

- {T 2 ~ a (0; T 2 (0 + ) + T 2 (0_)) ; + + 0_) . (A5) 

For the second iterate, we use 

T+'(l-ax;T 2 (T 2 (0 + ))+T 2 (T 2 (0_))) = I + 0(a), 
T+(l-ax;T 2 (T 2 (0 + ))+T 2 (T 2 (0_))) = 1 - ax' + 0(a 2 ), 
T 2a / (l-ax';T 2 (0 + )+T 2 (0_)) = -\ + 0(a), 



with 



and, finally, 



T 2 "(l -ax';T 2 (0 + ) +T 2 (0_)) = -1 + ax" + 0(a 2 ) (A6) 



x' = ^ + ^(T 2 (T 2 (0 + ))+T 2 (r 2 (0_))-2), 

= j + ^(T 2 (0 + )+T 2 (0_)-2) (AT) 



2^2a(l - a) v /x/16 - r 2 2 (0 + ) - r 2 (0_) 



The integrand, therefore, is given by 

d4> 



- (T+ (T 2 - (0; T 2 (T 2 (0 + )) + T 2 (T 2 (0_))) ; T 2 (0 + ) + T 2 (0_)) ; + + 0_) . (A9) 
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B The model function for p — > oo 

Because the second part of the radicand in Eq. (1221) contains a A;-fold Chebychev map, one first 
has to unfold this map. Using 

4>± = sin (^t±) , Po{4>±)d<P± = 2 d *±> ( B1 ) 



one obtains 

t: tn) = 



/^(At;to) = 7 f\dt + dt_x (B2) 



x cos |^7r 2 (A* - 2-Pad( fc )(t ;t+,t-)J -2arf (-cos(2 fc - 1 7rt + ),-cos(2 fe - 1 7rt_))j . 

After the substitution t' ± = 2 k t± the result reads 

f^\Akt ) = ^J^dt' + dt'_x (B3) 

9 

7T , N ,7T 



> cos j \/ | A/ - 2-''(/(/ 1 ''. ! (f ll :2-'f.2-¥j ) - 2nif ( - cos(-t' + ), - cos(-t'_; 

The huge integration range [— 2 k , 2 k ] for each of cf)' ± can be divided up into the standard intervals 
[—1, +1] by using 

4 = 4 + 2 fc -2n±-l G [-1,+1], (B4) 

where n± = 0, 1, . . . , 2 k — 1 counts the intervals. Because of t± = 2~ k (t± + 2n_(- + 1) — 1, the 

(two-fold linearized) draft function S k \to;t + ,t-) is nearly constant in each of the intervals 
[— 1,+1]. In defining 

d n+n _{t ) ■= lim rf (/c) (to;2- fe (2n + + 1) - l,2- fe (2n_ + 1) - 1) (B5) 

one has to deal only with the blunt. Using 

- cos(|4) = - cos - f (2 fc - 1 - 2n ± ] 

= - cos f |t" - 2 k ~\ + ^ + 7rn±J = sin Qt± + vm d 
= (-l)" ± sin(^) (B6) 
and removing the sign by replacing t± — > —t±, for the limit k — > oo one obtains 

l 2 "~ l r 1 



x 



n±=0' 



x cos I Jn 2 (At-2-Pad n+n _(to)Y -2ar? ( sin(^), sin( Jfr 



1 /" 

}™^7fc H / Po(0+)#+Po(0-)#- x 

n±=0 

^vr 2 (At - 2-%rf n+n _(t )) 2 - 2ar 2 °° (0 + , . (B7) 



n±=0 

x cos 1 



29 




In the final step, <p± = sin(7r£±/2) is used. The integrand does depend on ± only via r|°(0 + , </>_) 
while for the first part one can use the abbreviation 

c = it (At - 2~*ai n+n _ (t )) • (B8) 

In Fig. [TS] the function 

g<*\a, c) := J p o (0+)#+Po(0-)#- cos (y/& - 2ar p 2 (0+, (B9) 

is shown as a function of c for different values of a. For the order p of the blunt function 
7*2 (0+, 0~) a value of p = 3 is enough to provide sufficient precision. In analysing the position 
of the drafted minima in dependence on a, one obtain^ 



c(af - c(0) 2 = O 6(0) = ^c(a)2 + ^, (BIO) 
leading to 

<7«(o,c) = cos ^ + y I- (BH) 
This function can be inserted into Eq. ( IB 71) to give 



/»>(Af;to) = hm 1 cos f^r* (Af - 2-i>ad n+n _(t )) + y ] . (B12) 

In the limit k — > oo, the remaining summations over n± can be reunited to an integration over 
4>±, leading to the final result in Eq. (f23|) . For practical reasons, however, k is kept finite in the 
main text. 



*The value —8/3 is not rigorously derived but follows from numerical inspection. 
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C Combinatorics for the exponential draft function 



The (positive) coefficients of the exponent of the summation ordered exponential draft function 
can be calculated by using combinatorial considerations. In this appendix we demonstrate the 
procedure for single and multiple trees in the summation ordered expression and show how to 
derive the coefficients of the non-ordered expression. 



C.l Coefficients for single trees 

According to graph theory, single binary trees represent unique relations between powers of 2. 
The sum dfdf +1 di +s for instance, represented by the single tree 




results from the unique relation 1 + 1 + 2 + 2 + 2 — 8 = 0. In general, from combinatorics one 
can see that for the expression 

k— ri+l n—2 n—2 

E 11^ E"i = " ( C1 ) 

1=0 i=0 i=0 

there are n\/(no\ni\ ■ ■ ■ n n ^-) discriminable permutations without taking into account the global 
sign (which gives a factor 2!). Together with a factor 1/n! for the exponential series in Eq. ( I33p 
and a factor l/2 n from the trigonometric functions T 2 !(sin(7rt + /2)) after their expansion in 
exponential functions, one obtains 



2\n\ 



n\2 n \n !ni! • • -n n _ 2 !/ 2 n 1 n !ni! • • • n n __ 2 ! 



(C2) 



C.2 Coefficients for multiple trees 

While the coefficients of the single tree contributions are the same for the summation ordered as 
for the non-ordered exponential draft function, the coefficients of the multiple tree contributions 
are different in both cases. However, the coefficients of the summation ordered contributions 
can be calculated by combinatorial means. This is shown for the contributions of orders n = 4 
to 6 in the following. Multiple binary trees correspond to relations between powers of 2 which 
can be split into two or more single relations separated by a vertical bar. 

• £df: +1-1|+1-1 = gives 4!/2!/2! permutation if +1 and -1 are taken to be 
two different elements. A global sign for the first or the second part +1 — 1 = need 
not be taken into account because the permutations already lead to the sign-changed 
contributions (sign non-sensitive). Together with the global factor 1/41/2 4 one obtains 

1 / 4! \ 1 1 
412 1 \2!2! J ~ 2 4 2!2! ~ 64' 
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E^M+i : +1 + 1 — 2| + 1 — 1 = gives 5!/3!l!l! permutations where now the global 
sign for the first part has to be taken into account. One obtains two (formally equal) 
contributions 

1/5! 5! \ 2 1 

5!2~5 ^3!1!1! + l!3!l! J ~ 2^3! ~~ 96' 

J2dfdf +1 : +1 + 1 — 2| + 2 — 2 = gives 5!/2!l!2! permutations where the same holds as 
in the previous case. One obtains 

1/5! 5! \ 2 1 

5!25 V2!l!2! + 2!2!1! J ~ 2 5 2!2! ~ 128' 

Edfdf +1 di +2 : +1 + l-2| + 2 + 2 -4 = gives 6!/2!l!2!l! permutations. In this case one 
has to take into account the global signs of both parts. In addition the numer of elements 
changes if the global signs change. One obtains 

1 / 6! 6! 6! 6! \ J_ 

6!26 l v 2!l!2!l! + 2!3!1! + 2!1!2!1! + 2!3!1! J ~ 96' 

E dfd l+1 d l+2 : +l-l| + l + l + 2- 4 = gives 

1/6! 6! \ 1 

6!26 \m\m\ + 1!3!1!1! J ~ 192' 

E dfd l+1 df +2 : +l + l + 2-4|+4-4 = gives 

1/6! 6! \ 1 

6!26 l v 2!l!l!2! + 2!l!2!l! ) ~ 128' 

E dfdf +1 : Besides the sign sensitive possibility +1 + 1 — 2| + 1 + 1 — 2 = 0, in principle there 
is also a non-sensitive possibility with three parts, +1 — 1| + 1 — 1| + 2 — 1 = 0. However, 
this possibility is the same as the sign-changed possibility +1 + 1 — 2| — 1 — 1 + 2 = 0. 
Therefore, only the first has to be taken into account, resulting in 

1/6! 6! 6! 6! \ 1 

6!26 ^412! + 2!2!1!1! + 2!2!l!l! + 4!2! ) ~ 1536 ' 

E^f : +1 — 1| + 1 — 1| + 1 — 1 = contains three parts not sensitive to the global signs. 
Therefore, one obtains 

— (-) = — 

6!2 6 \3!3! J 2304 
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C.3 Summation ordered and non-ordered coefficients 

The relation between the coefficients of the summation ordered and non-ordered exponential 
draft function can be found by using 

£<* 2 £4 = + 



3 

Z+2' 



E dfdi + idi +2 E d 2 , — £ dfdi + id l+2 £ c^, + £ dfd[ + idi +2 + £ d 2 df +1 d l+2 + £ d 2 di +1 d] 
£d / 2 <iz + i£<if,d Z / +1 = 

= £ <i^ + i £ d 2 v d v+ i + £ d 2 df +1 d l+2 + £ <^H 2 +i + £ d 2 f df, +1 d[' +2 = 
i 

= £ d l di + i £ d v dy + i + £ di d l+1 + 2 £ <^ d l+1 d i+2 , 

£4£4 = £#£*# + £<£ 

£4£4£4 = (£4£4)£4 + £4£4 = 
= E^ 2 E4E4 + E^ 2 E4 + E^E4 + E4E4 + E^ = 
= E4E4E4' + 3E4E4 + E4- (cs) 

Using, for instance, for the contribution of order n = 6 in Eq. ( I34p a general ansatz and applying 
Eqs. (lC3p . one obtains 

^i E d^di + idi +2 di +3 di +4: + A 2 £ d 2 df +1 d[ +3 + A 3 £ dfdi +2 di 3 

+ At £ dfdf +1 d l+2 + A 5 £ dfd l+1 d h + A 6 £ d 2 di +1 d 2 l+2 + A 7 £ df d 2 +1 

+ A 8 £ d 2 di +1 d i+2 £ t$ + A 9 £ d 2 G^+i E + ^io E ^ 

+ ^iiE4E4 + ^ 2 E4E4E4 = 

= ^i E dfdi + idi +2 di +3 di + 4, + v4 2 E dfdf +1 d[ +3 + A 3 £ dfd[ +2 di s 
+ (A 4 + A 8 + 2A 9 ) £ d 2 df +1 d t+2 + (A 5 + A 8 ) £ df d /+1 d, 2 
+ (A 6 + A 8 ) £ d 2 d i+1 d 2 +2 + A £ dfd 2 +1 + A 8 £ d 2 d l+l di +2 £ 4 
+ A 9 £ d 2 d l+1 E d 2 v d u+l + (Am + A n + A X2 ) £ df 

+ (An + 3A 12 ) £ di £ 4 + ^12 E d] £ 4 £ 4- (C4) 
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Comparing with Eq. (j4"0|) one obtains 

1 1 1 a 1 1 
Ai = — , A 2 = , A 3 = , A 8 = — , A g = . 

64' 384' 768 64' 128' 

furthermore 

1,11 1 

An = 3v4i2 = = , 

256 256 128 256 

and finally the coefficients for . . . 

^2dfdf +1 d l+2 : At 
53 dfd l+1 di +2 : A 5 
^2dfdi + idf +2 : A 6 



1 

96 ~ 


-A 8 - 


-2A 9 = 


1 

" 96 ~ 


1 

" 64 


1 

- 2— 
12, 


1 

192 


- A 8 


1 

~ 192 


1 

' 64 




1 

96' 


1 


-A 8 


1 


1 




1 


128 


" 128 


' 64 




128' 



1536 1536 128 1536 

V df : A w = — A u - A 12 = — h — 

^ 1 2304 2304 256 3 
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